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Abstract 

Unboundedness of solutions is shown to occur in a one-dimensional quasilinear parabolic- 
parabolic chemotaxis system for any initial mass. Our result is also independent of the relation 
between the speeds of the diffusion of cells and chemoattractant. The proof is achieved by 
contradiction. It uses a virial type inequality together with a boundedness from below of the 
Liapunov functional associated to the system. Moreover, the same method is applied also in 
higher dimensions and an unboundedness result is shown in the case of critical quasilinear 
fully parabolic Keller-Segel system for mass large enough in dimensions n > 3. 
Key words: chemotaxis, infinite-time blowup. 
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1 Introduction 

This work deals with nonnegative solution couples (u, v) of the parabolic-parabolic Keller-Segel 
system 

u t = (a(u)u x ) x - (uv x ) x , x € (0, 1), t > 0, 

rv t = v xx —v + u, x e (0, 1), t > 0, 

u»(0)=u x (l)=« !B (0) = i; as (l)=0, t>0, 

u(x, 0) = uq(x), v(x, 0) = vo(x), £ G (0, 1), 

where t > and the initial data are supposed to satisfy uq, vo € W 1,p (0, 1), p > 1 such that uo > 
and vq > in (0, 1). Moreover, we assume that a £ C 2 ([0, oo)) is a positive function. 
Systems of this kind were introduced in [T5] to describe the motion of cells which are diffusing 
and moving towards the gradient of a substance called chemoattractant, the latter being produced 
by cells themselves. The original model consisted of several equations including more unknowns. 
It was simplified in [19] where the form of the equation similar to (jl.l[) was introduced. The 
essentiality of the nonlinear diffusion was emphasized in |llj where the so-called volume- filling effect 
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was studied. The important issue in studies concerning (jl.ip was a chemotactic collapse of cells 
interpreted as a finite-time blowup of the component u. In 3 , based on the numerical simulations, 
the authors suggested that in a two-dimensional (most popular and biologically most relevant) 
scmilinear version, there is a threshold of mass such that a solution starting from the initial data 
with a mass smaller than a threshold exists globally, thus avoiding the chemotactic collapse, while 
initial conditions having mass larger than a threshold lead to a chemotactic collapse. We emphasize 
at this point that the rigorous studies that followed were related to further simplifications of the 
Keller-Segel model. Actually, until recently, almost all the results concerning finite-time blowup 
were achieved for a parabolic-elliptic simplification (r = in First, in P3] a conjecture 

from [3] was shown in the two-dimensional semilinear parabolic-elliptic simplification. In j!6) the 
exact threshold in the case of radially symmetric data was calculated as 8tt in the parabolic-elliptic 
2d case. Next, the critical mass phenomenon was also proved in the case of a two-dimensional 
problem for solutions without symmetries (see [17], [18]). Moreover, it has been shown that in 
higher dimensions a finite-time blowup of solutions to the semilinear version of (|1.1|) can occur 
independently of the initial mass provided that the initial data are concentrated enough [16]. 
Finally, during last years, the quasilincar system was studied. Critical nonlinearities (and the 
lack of it in dimension 1) and critical mass phenomenons were found in all the dimensions, see 

However, all those results are available only for a parabolic-elliptic simplification of (jl.ip . In the 
case of the original fully parabolic version the investigation of chemotactic collapse turned out 
to be a much more challenging issue. Until 2010 the only available result stating the occurrence 
of finite-time blowup was the one in [12] . where examples of some particular solutions to the 
semilinear version of (jl.ip in dimension n = 2 blowing up in a finite-time are shown. Moreover, 
see [13], [21], it was known that solutions to a fully parabolic problem can be unbounded. However 
it was not stated if blowup occurs in finite or infinite time. The same approach was used in 
[22] to show unboundedness of solutions to a fully parabolic Keller-Segel system in the case of 
supercritical quasilincar systems in dimensions n > 2. Next, in [5] the collapse of solutions to the 
one-dimensional Keller-Segel system with appropriately weak diffusion of cells and sufficiently fast 
diffusion of chemoattractant is shown provided initial mass is large enough. The breakthrough has 
been made recently in [25] , Introducing a new method M. Winkler shows there that in dimensions 
n > 3 generic solutions to the semilinear version of (jl.ip blow up in finite time independently of 
the size of initial mass. His method was generalized to the quasilinear case in [5], [5], where the 
finite-time blowup of solutions to a fully parabolic supercritical Keller-Segel system was shown for 
any initial mass in any dimension n > 2. The problem for critical nonlinearities is still open. 
This way, the occurence of finite-time singularities of solutions to a fully parabolic quasilinear 
Keller-Segel system in higher dimensions is already known for supercritical diffusions independently 
of the size of initial mass. However, in the one-dimensional case, the only known result is the one 
from [5] , where the finite-time blowup is shown only for large mass initial data and moreover in a 
neighbourhood of a parabolic-elliptic situation, meaning r being a very small number. The purpose 
of the present paper is to introduce a result concerning the supercritical case of (jl.ip with any r > 
and an arbitrary mass. We are going to give a result stating the unboundedness of solutions. It does 
not imply a finite-time blowup. We do not exclude a situation that a solution becomes unbounded 
when time goes to infinity. We know that this kind of situations might happen (see [5] Theorem 
1.6] or [9] Theorem 1.4]) in chemotaxis. Our result is similar to the higher-dimensional results in 
[21], [13] or [22], however the method of a proof which we follow is totally different. The method 
of [21] , [13] , [22] is not applicable in the one-dimensional case. Indeed, unboundedness of solutions 
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of a higher dimensional Keller-Segel system in the above-mentioned articles is a consequence of a 
contradiction between boundedness of solutions which immediately leads to convergence with time 
(at least on subsequences) to steady states and the fact that one can construct initial conditions 
such that the value of the Liapunov functional on them is smaller than the minimum of this 
Liapunov functional over the set of steady states. This approach fails in the one-dimensional case 
due to the fact that the Id Liapunov functional is bounded from below by its minimum over the 
steady states. We develop an alternative way of proving unboundedness of solutions in that case. 
We use an extension of a moment method (see [2J, Qj5] or [5]). 

The same method can be applied also in higher dimensions. However, in supercritical cases, 
unboundedness results, as well as much stronger finite-time blowup results, are already known, 
see [T3], [3T], [52], [13], [5]. But so far, there were no results (neither those concerning finite nor 
infinite-time blowup) in the critical case in dimensions n > 3. We fill in this gap in the last section. 



2 Preliminaries 

First, notice that thanks to the application of Amann's theory ([T]) we have a local well-posedncss 
result for the regular nonnegative solutions, once we start from the regular nonnegative initial 
data, see ([!]). Integrating by parts we see that those solutions also satisfy 

u(t,x)dx = j uo(x)dx, / v(t,x)dx — / uo(x)dx + e~~ I / vo(x)dx — / uo(x)dx 
Jo Jo Jo \Jo Jo 

(2.1) 

for any t > 0. Moreover, if the solution ceases to exist for times larger than to, then immediately 
L°° norm of u blows up at a time to- 

Next we recall the existence and properties of a Liapunov functional associated to (jl.ip . 

L{u{t),v{t))+r f \\d t v{s)\\l ds+ f [ (tt( " K " mil)2 =L(u ,v ) for t>0, (2.2) 



^0 Jo 



where 



1 ,„ ,o 1 



II 



L(u, v) := / ( *(«) - uv + ^ |a^| 2 + - M 2 ) dx , (2.3) 



where $ G C ((0, co)) is defined in a following way 

$(1) = $'(1) := and $"(r) := ^ for r>0. (2.4) 

r 

We take an advantage of the one-dimensional setting and state the boundedness from below of the 
Liapunov functional, compare [5] Lemma 5] 

L(u(t),v(t)) > -C for t>0. (2.5) 

The main theorem of this part reads 

Theorem 2.1 Assume a 6 -^(O, oo). Then for any r > and any initial mass m > there exists 
an initial data (uo(x),vq(x)) such that J uo(x)dx = m and a component u(t) of a solution to (jl.ip 
starting from (uq,vq) is unbounded in an L°° norm. 
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Remark 2.2 Our theorem does not say whether the blowup occurs in finite or infinite time. 

In the sequel we will need the following observation which is a consequence of classical regularity 
estimates for parabolic equations, however in our case we can infer it easily, keeping the proof 
self-consistent. 

Proposition 2.3 Assume \\u(t, OIIl^o i) < @ f or t > 0, then there exists a constant C2 > such 
that \\v t (t, -)|| 2 < c 2 fort > 0. 

Proof. First we differentiate the second equation in (jTTTJ) in t and next multiply it by Vt and 
integrate. We arrive at 

Tj-jT / v 2 dx + I (v tx ) 2 dx+ [ v\dx = [ u t v t dx. 
£ at J J J J 

By the first equation in we see that 

u t v t dx = / (a(u)u x - uv x )v tx dx < - [ (a(u)u x - uv x ) 2 dx + i / (v tx ) 2 dx. 



t d f 1 2 1 f 1 2 f 1 2 , , IMIl~(o,i) f 1 {a(u)u x ~ uv x ) 



Hence 



Next we observe that due to (|2.2|) and (|2.5|1 we have 
1 t f 1 (a(u)ux - uv x ) 2 



v 2 dx + -J (v tx ydx+ I v 2 dx< / .. — dx - 



n 



\2 



^0 



dxcii < oo for any i > 0. 



Thus, integrating (|2.6p with respect to time, keeping in mind that \\u(t, Olliooro i) ^ s assumed to 
be bounded, we arrive at the desired estimate of v t . □ 



3 A proof of the main result 

This section is devoted to the proof of Theorcm l2.ll We are going to use a method of a generalized 
second moment as introduced in [5]. However, our aim is weaker than in a mentioned paper. 
Thanks to the estimate on v t coming from Proposition 12. 31 by applying some ideas from [5] , we 
arrive at the estimate which in turn allows us to infer unboundcdness result for arbitrary initial 
mass and arbitrary r > 0. This result is of special interest due to the fact that among supercritical 
Keller-Segel systems, the one-dimensional one was by now the only case where not only finite-time 
blowup, but even unboundedness of solutions was not known. In dimensions n > 2 it is already 
known in supercritical cases that for any mass there exist solutions blowing up in finite time, see 
[5], [5]. However, the method used in those cases, being an extension of the method of Winkler, 
[23] , does not seem to be applicable in the one-dimensional case due to the bound of the Liapunov 
functional from below, see (|2.5[) . For the same reason the usual method of showing unboundcdness 
of solutions to the higher-dimensional Kcllcr-Scgel system does not work in the Id setting. This 
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method, see [13] , [S] , H2 , requires unboundedness af a Liapunov functional in order to enable the 
user to choose initial data such that the value of the Liapunov functional on them is smaller than 
the innmum of the values of the Liapunov functional over the steady states. In a onc-dimensional 
setting the minimal value of the Liapunov functional seems to be attained at the steady states, 
excluding this way the above method. 
Our proof is of a reduction to absurd nature. 

Proof. Let q > 2 and define the generalized moment in a following way 



M q (t) 



1 Jo 



\U(t,x)\ q dx, 



where 



pX pX 

U(t,x) := / u(t, z)dz for x G [0, 1] and V(t, x) := / v(t, z)dz for x € [0, 1]. 
Jo Jo 

Assume m := J Q uo(x)dx > 0, 



0< / vo(x)dx — m < 



2(9 + 1)' 



Moreover assume sup t>0 \\u(t, -)lli,=°(o l) < 00 • By (QUI) we have 



j t M q (t) = I Ui-^Aiu))* - I U^uirVt -U + V), 



(3.1) 



(3.2) 



where A(u) := — J°° a(s)ds. Integrating by parts we have 
d 



dt 



M q (t) = -(q-l) U g - 2 uA{u) + (U q - 1 A(u)) + U q - X u{p -V -rV t ) 



(g-1) / U q - 2 (-uA(u)) + m q - 1 A(u(t,l)) + 



1 



9 + lio 



(U q+1 ) x 



q Jo 



{U q ) x {V + rV t ). (3.3) 



Due to (f2~Tj) and ([33]) we have 



9 Jo 



1 Jo 



{U q ) x V = - U q v- -(m q V(t,l)) < - / U q v 



1 Jo 



- m i +1 



which in turn, together with (|3.2I) and an observation A(u(t, 1)) < yields 



-M q (t)<(q-l) ^ t^- 2 (-«A(ti))+||«|| ioo(0tl) Af^Ct)--^-^ 
Since 



,9+1. 



U q v t ) -m q V t (l) 
(3.4) 



= -e" 

T 



i>ocfcE — ml, 



■5 



d 
dt 



Next, by Holder's inequality 



,<?+! 



in view of (|3.ip we have 

M q {t) < j U«-\-uA(u))+\\v\\ Lxm MS)-^^+^ U«v t ) + 

<(q-l) [ Ui- 2 (-uA{u)) + \\v\\ Laom M q {t) 



2q{q+l) 



dt 



M,(t)<(g-1) / U"- 2 (-uA(u)) + \\v\\ L ^ {0A) MS)- 



1 



m q 

q 

U"vt 



Vodx — rn 



,9/ 2 



2qiq + r) m^ + —M ?M2 .(3. 5) 

We are now in a position to introduce a concave function B such that 

< -rA(r) < B(r), r G (0, oo) , 



lim — ^ = 0. 



(3.6) 
(3.7) 



This function was introduced in [6| and in [7j Lemma 3.1] it was shown that such a function can 
be associated to any a satisfying the assumptions of Theorem 12.11 Its advantage is that we can 
proceed for any m > 0. The following estimates, using in a tricky way Jensen's inequality for 
several probabilistic measures, are taken from (6j Theorem 10], we provide them for the sake of 
completeness. 



/ U q - 2 B{u)dx 
Jo 



B{u)dx) / (U q ) — — 



s=2 B(u)dx 



Jo B{u)dx 



< 



B{u)dx 



l uq _ B(u)dx 



o Jq B(u)dx t 



B(u)dx (qM q )~ 



< (qM q )^ 







\ B ( f 1 uU"dx\l 












I Wo lM q )\ 



1 ft( \ uqdx 

B(u) 

V ' qM q 

< (qM q )^ [B{m)f/ q 



B 



,9+1 



q(q + l)M q 



The Jensen inequality above was applied to j^g|^J » ^ x an< ^ U qM X • 

In view of (|3.8[) . Proposition 12.31 and the bound on |H| ioO ( i)i <|3-5p reads 



(3.8) 



-M q (t)<(q-l)(qM q )^[B(m)] 



Next ((321) reads 



1 uU q dx\ 



o 9M, 



-J 



,0+1 



7?/ 2 



-M (t) < A ro (M (t)), 



ClM ^-2^Ti) + ^ M - 

(3.9) 
(3.10) 
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where 



9+1 

A(r) := Cl r + (q - l)[B(m)] 2/q (—— 

\q+l 

for r > and /3(r) = 

We notice that according to (|3.7|) for any to > A m (0) < 0. Due to continuity of A, there exists 
6 > such that A(r) < for r e [0,0]. Thus, being allowed to choose initial data such that 
M q (0) < 6 we arrive at the claim that M q must touch at a finite time. But this is contradictory 
to the nonnegativity of u and a mass conservation (|2.1[) , what means that a solution ceases to exist 
at latest at the time of an extinction of M q . This in turn means that u is not bounded in a 
contradiction. 

□ 



,9+1 



q(q + l)r 



C2TTO 



<//- 



'1 



1/2 



-7-2 - 



2q(q+l) 



4 An application of the method in higher dimensions 



This section is devoted to proving unboundedness of radially symmetric solutions to the fully 
parabolic quasilinear critical Keller-Segel system in dimensions n > 3. We achieve this issue 
extending the method of previous section to higher dimensions. Again, our result does not say 
whether blowup of solutions occurs in finite or infinite time. 

As is well known (see for example |10| ) radial symmetry is preserved by the Keller-Segel model in 
a ball and equations of evolution in that case read 



Ut = r^r(r" 1 a(u)u r ) r - -^=r(r n 1 uv r ) r , 



TV""" 
1 f„n—\. 



r e (0,1), t > 0, 
v t = pr=r(r n - l «r)r +u-v, re (0, 1), t > 0, 
u r = v r = 0, re {0, 1}, t > 0, 
u(r, 0) = uo(r), v(r, 0) = ^o(r) r £ (0, 1). 



(4.1) 



Like in the case of u, u denote the density of cells and chemoattractant. Local-in-time solutions 
exist and whenever they cease to exist for all times we know that \\u(t, ■)IIl°°(b(o i)) S oes to infinity 
for some t < oo. Moreover both uo,vq are assumed to be nonnegativc, hence for t > both 
u(t, -),v(t, •) are positive due to a comparison principle. The function a is C 2 regular and positive, 
it stands for nonlinear diffusion. The radius of a ball is 1, and the radial variable is denoted by 
r : < r < 1. Next, 



Jb(i,o) x "> dx = /u(i,o) Mx)dx, 
/b(i,o) = is(i,o) «o(a;)da; H 

hold, where B(0, 1) is a ball centred at 0. 
Like in an one-dimensional case, we have 



(Jb(i,o) M x ) dx - Jb(i, ) "0(^)^2; 



(4.2) 



|9tu(s)||2 rfs 



|a(u)Vu — uVd | 2 



(0,1) 



i(«o,«o) for f>0, (4.3) 
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with L defined exactly as in (|2.3[) with integration over the interval replaced by the one over a ball. 
Moreover, we have a simple proposition 

Proposition 4.1 Assume \\u(t, ')IIl~(b(o i)) < @ f or t > 0, then there exists C\ > such that 
L > -Ci. 

This is implied by the fact that J B r ^ uv is bounded independently of time. Next, analogously to 
Proposition 12 . 31 we have 

Proposition 4.2 Assume \\u(t, -)IIl°o(_b(o i)) < ^ f or ^ > ®> then \\vt(t, -)|| 2 < C\ for t > 0. 
We are now in a position to formulate the main theorem of this section. 

Theorem 4.3 Let n > 3. Assume a(u) = (1 + TTien there exists such value, say m*, 

that for any M > to* there exists (uo,Vo) with mass J B r j\ Uo(x)dx = M\B(0, 1)|, smc/i i/iai a 
component u(t) of a solution of (|4.1[) starting from u$ is unbounded in L°° norm. 

Proof. Like in section [3] first we assume that u is bounded in L°°. In what follows we show that 
this assumption leads to an absurd. Let us follow the ideas which appeared in pQ. We introduce 



U(t,r) :-- 



u(t,x)dx and M 2 (t) 



n\B(0,l)\J mr) 
for (t,r) £ [0,T max ) x [0,1]. Next, we notice that 



M 



U(t,0) = U(t,l) = 0, 



(4.4) 



Moreover for V defined analogously to U we have 



M 



V(t, 1) = — + e-* V(0, 1) - 



M 



(4.5) 



Integrating (|4.1[) we get 

C/ 4 = r"- 1 A(w) r - r™" 1 ^, F t = r™~ V - V + E7, 
where A is a primitive of a such that A(0) = 0. Hence 

dM 2 f 1 (M \\ „,„ , M (M \ 

-df = I U " u ) r " + (~ ^) ^ + yc/ " + Kt 

Consequently, 

rfM 2 



?7r 



dr. 



2(n - iy n ' A \ —-U)- 



M 



^(n-l) 



U, 



At \1 1 f A/ 

d V n 









f 


} 
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Thus, due to the fact that integration by parts gives 

2 















/o 











we infer 



dM 2 
dt 



< 



L 



AI 



2(n - l)r 2 "" 3 17 - r 2 ("-^C/ r 



1 / M 



A(u)dr- - I — ) + IMI L ~ (B(0 ,i)) 



M 2 (t)+ 
(4.6) 



M 



(Ma(t))' / <r n - x dr 



.-/() 



M 



At ^ J 1 f M 
i4(u)dr - - — 
6 V n 



+ IMIl~(B(0,1)) M 2 (t) + n3/2|Jj(M)|1/2 (M 2 (f))^ IMII W 1)) 

Keeping in mind that a(u) = (1 + , we see that 



A(u) < 



.2-- 



2- ^ 



Hence we are now in a position to use computations from [3] in order to estimate a term 



AI 



2(n-l)r 2 ™- 3 — -U) -r 2( - n ~Vu, 



A(u)dr 



For the sake of completeness we provide them with all details. First we notice that 



(4.7) 



2(n -1)1— -U 



A(u) < max{2(n -1)[—-U)- r n u, 0} 



1 



2- H 



Next, 



[ — -U 



1 



max{2(n — 1) 
1 



U - r n u, 0} 



2- 2 

n 

in turn (14. 7[) is estimated by 
/■l 



M 



2(?i — 1) ( — — C7 



1 



V 2 ~ 

2— a 



2 - u 
2(n- 1) / M 



-it 2 " + u | < 



r" n+ ^ + 2(n-l) [ — - ET ) u, 



17 



1 \ u < 



1 



'0 



2- 2 



2(n- 1) J7 



A/ 



C/ r dr + 2(n - 1) / r 



M 



U U r dr 
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A term 2(n - 1) / r™~ 2 (M _ u) U r dr can be estimated further using the fact that the function 
r — > is concave. We use the Jensen inequality applied to a measure — U) U r dr, integrate 

by parts and obtain 

2(n-l)J r n - 2 (j£- u J U r dr< (n-l)(2n) 1 "« ffl"if a (t)H. (4.9) 
We conclude that in view of (|4.8p and (|4.9I) , (|4.7j) is estimated by 

(2(n-l)) 2 ~t /M\ 3 "» , „,,„,, >/M^„ n , 2 
J +(n-l)(2n) 1 -- f — J M 2 (ty-«, 



(2-S)(3"f) 

hence (|4.6p yields 

"\2-# X 71 ;T X 3 "ft 1 /H\ 3 



dM 2 (i) (2(?i-l)) 2 -- [MY " 1 /M\ , w x1 ! /M\ " , 2 



M 



(B(0,1)) M 2(*) + n 3/2| B ( 0il )|l/ 2 ( M2 W) 3 II U *II1 = (B(0,1)) 



Now, due to an assumption that L°° norm of u is bounded and a local existence result, we infer 
that solution exists for any time t > 0. However, as we see in a moment, for initial mass greater 
than to* we can find such initial data (uq,Vq), that this is in contradiction with (|4.10[) . Indeed, 
according to Proposition 14.21 11^11^2(3(0 i)) ^ s bounded, moreover v is also bounded in L°° norm 
due to the standard parabolic regularity estimates related to the second equation in (|4.ip . Hence 
for M large enough that 

(2(n-l)) 2 ~* (M\^ 1 /AT " 



(2-f)(3-|) 6 

and choosing initial data uq concentrated enough that M2(uo) is close to 0, we see that the right- 
hand side of (|4.10[) is negative at the initial time. But this means that M2{t) is a decreasing 
function with a negative time derivative bounded from above, so Mz(t) hits at a time to < oo. 
Due to nonnegativity of u and mass conservation, this is in contradiction with global existence of 
a solution. 

□ 



Remark 4.4 Notice, that analyzing a proof of Theorem \4-3\ we conclude that once we pick up a 
proper Uq then whatever is a Vq we start from, anyway our solution must become unbounded. This 
is unlike in a one- dimensional setting of Theorem \2.1\ where vq had to be chosen carefully so that 
\3.1\) is satisfied. 
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